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Abstract The present paper explores the synchroniza- 
tion scenario of hyperchaotic time-delayed electronic os- 
cillators coupled indirectly via a common environment. 
We show that depending upon the coupling parameters 
a hyperchaotic time-delayed system can show in-phase 
or complete synchronization, and also inverse-phase or 
anti-synchronization. This paper reports the first exper- 
imental confirmation of synchronization of hyperchaos 
in time-delayed electronic oscillators coupled indirectly 
through a common environment. We confirm the oc- 
currence of in-phase and inverse-phase synchronization 
phenomena in the coupled system through the dynam- 
ical measures like generalized autocorrelation function, 
correlation of probability of recurrence, and the concept 
of localized sets computed directly from the experimen- 
tal time-series data. We also present a linear stability 
analysis of the coupled system. The experimental and 
analytical results are further supported by the detailed 
numerical analysis of the coupled system. Apart from 
the above mentioned measures, we numerically com- 
pute another quantitative measure, namely, Lyapunov 
exponent spectrum of the coupled system that confirms 
the transition from the in-phase (inverse-phase) syn- 
chronized state to the complete (anti-) synchronized 
state with the increasing coupling strength. 



Keywords Delay dynamical system • hyperchaos • 
chaos synchronization • environmental coupling • time 
delay electronic circuit 



T. Banerjee* • D. Biswas 

Department of Physics, The University of Burdwan, Burdwan 
713 104, West Bengal, India. 

* E-mail: tanbanrs@yahoo.co.in, tanban.buphys@gmail.com 



1 Introduction 

For more than two decades, synchronization of chaos 
and hyperchaos have been an active field of research in 
various fields, including physics, biology, mathematics, 
engineering, etc. After the seminal paper of Pecora and 
Carrol pQ on chaos synchronization, several synchro- 
nization schemes and processes have been observed and 
identified; example inludes, complete synchronization, 
generalized synchronization, phase synchronization, lag 
synchronization, anticipatory synchronization , impul- 
sive synchronization, etc. (for detailed discussion re- 
fer to two excellent review papers [2j[3] and references 
therein). 

After the initial emphasis on the synchronization 
of low dimensional chaotic systems (i.e., chaotic sys- 
tems having a single positive Lyapunov exponent), re- 
searchers soon attracted towards the study of synchro- 
nization in high-dimensional systems like complex net- 
work [4], and nonlinear delay dynamical system (DDS). 
Particularly, studies on the synchronization scenario of 
delay dynamical systems are challenging, both from 
theoretical and experimental point of view, owing to 
their infinite dimensionality with a large number of 
positive Lyapunov exponents. Further, in real world 
most of the time we have to encounter with time-delay; 
few examples are, blood production in patients with 
leukemia (Mackey-Glass model) [5], dynamics of optical 
systems (e.g. Ikeda system) [6 [7 , population dynam- 
ics [8], physiological model [9], El Niho/southern oscil- 
lation (ENSO) [10], neural networks [11 , control sys- 
tem p~2|n~3| fT4] . etc. The first study on synchronization 
of chaos in time-delayed system has been reported by 
Pyragas [15 j. Later, several genres of synchronizations 
in time-delay systems have been reported; few examples 
are (without claiming to be complete): Lag and antici- 
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patory synchronization [16^ 7, 18, 19, 20 , complete and 
generalized synchronization [21,22,23 , etc. The phe- 
nomenon of phase synchronization in time-delayed sys- 
tems has first been reported by Senthilkumar et al [24] 
and later it was confirmed experimentally in Ref. [25] . 
Recently, global phase synchronization and zero-lag syn- 
chronization (ZLS) in time-delay systems have been re- 
ported in [26] and [27] , respectively, showing the ongo- 
ing interest in the field of synchronization of time-delay 
systems. In all the above mentioned works (except [27]), 
the coupling scheme is essentially the direct coupling, 
i.e., either unidirectional coupling or bi-directional cou- 
pling, where either of the two coupled systems or both 
the systems directly affect the dynamics of each other. 
Ref. [27] considers mutual and subsystem coupling con- 
figurations via dynamical relaying and reports the ex- 
perimental confirmation of ZLS in a system of three 
coupled piecewise linear time-delay circuits. 

An important coupling scheme has recently been 
reported, namely the environmental coupling , which 
is an indirect coupling [28] scheme. Here two or more 
systems interact indirectly with each other through a 
common environment. The environment makes the sys- 
tems synchronous without affecting their central dy- 
namical features. Environmental coupling is particu- 
larly important in biological systems, e.g., populations 
of cells in which oscillatory reactions are taking place 
interact with each other via chemicals that diffuse in 
the surrounding medium [28 . The environmental cou- 
pling scheme has been investigated numerically in de- 
tail by Reshmi et al [29] where the phase and complete 
(anti-) synchronization of environmentally coupled low- 
dimensional chaotic systems (Rossler and Lorenz sys- 
tems) are explored using Lyapunov exponent spectrum, 
average phase errors, and correlation functions. A linear 
stability analysis has also been reported in that paper. 
Later, the experimental study of the same has been 
reported in [30] . Amplitude death in environmentally 
coupled chaotic systems are reported in [31,32;. In this 
context, another important work has recently been re- 
ported by Sharma et al. [33] that considers a variant of 
environmental coupling to show (both numerically and 
experimentally) the occurrence of phase-flip bifurcation 
in periodic and chaotic systems. 

All the previous works on environmental coupling 
are restricted to the study of low- dimensional chaotic 
systems. The same for the time-delayed hyperchaotic 
system has not been reported yet; particularly exper- 
imental realization of environmentally coupled hyper- 
chaotic time-delayed system is challenging enough, and, 
to the best of our knowledge, yet to be reported. In this 
paper we study the synchronization states of two hyper- 
chaotic time-delay systems that are coupled indirectly 



via a common environment. For the present study we 
choose a hyperchaotic time-delay system recently pro- 
posed in [34[. The choice of this system is led by the 
fact that unlike other time-delay systems it poses three 
important features that are useful for the experimental 
design: first, it has a closed form mathematical function 
for the the nonlinearity (unlike piece- wise-linear nonlin- 
earity), second, it shows hyperchaos even for a small or 
moderate value of time-delay, and lastly, and most im- 
portantly, it can be realized with off-the-shelf electronic 
circuit components. Further, the hyperchaotic attractor 
of the system is inherently phase-incoherent, thus the 
present study can be extended to other general class of 
hyperchaotic systems. We show that depending upon 
the coupling parameters a hyperchaotic time-delayed 
system can show in-phase or complete synchronization, 
and also inverse-phase or anti-synchronization. For the 
first time we report the experimental studies of syn- 
chronization of hyperchaos in time-delayed electronic 
oscillators coupled indirectly through a common envi- 
ronment. We confirm the occurrence of phase synchro- 
nization in the coupled system through the dynamical 
measures like generalized autocorrelation function, cor- 
relation of probability of recurrence [35,36^, and the 
concept of localized sets [37 computed directly from 
the experimental time-series data. We also perform a 
linear stability analysis of the coupled system for the 
complete and anti-synchronized cases. The experimen- 
tal results are further supported by the detailed nu- 
merical simulations of the coupled system. Numerical 
computations are carried out to find out the Lyapunov 
exponent spectrum of the coupled system that confirms 
the transition from the in-phase (inverse-phase) syn- 
chronized state to the complete (anti-) synchronized 
state with the increasing coupling strength. Numerical 
recurrence analysis and concept of localized set are used 
to reconfirm the occurrence of phase sysnchronization 
in the coupled system. 

The paper is organized in the following manner: 
the next section describes the environmental coupling 
scheme for a time-delay system, and also describes the 
system under consideration and a summary of its main 
dynamical behavior. Section [3] describes the experimen- 
tal design of the coupled system. Experimental results, 
recurrence analysis and other dynamical measures are 
presented in this section. Sect[4]gives an account of lin- 
ear stability analysis of the synchronized state. Sect. [5] 
describes the numerical simulation results of the cou- 
pled system, Lyapunov exponent spectrum, recurrence 
analysis, concept of localized set, and two-parameter 
stable zone of synchronization have been reported in 
this section. Finally, Sect. [6] concludes the outcome and 
importance of the whole study. 
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Fig. 1 Nonlinearity with the function f(x T ) = —nx T + 
mtanh(lx T ) with "Nl": n 2.2, m = 1, I = 10; "N2": 
n = 3, m = 1, Z = 8, "N3": n 1.5, m 1.2, Z 8. Dashed 
curve is for n = 4,ra=l,Z = 3, which shows that for n > ml 
the nonlinearity does not show the two-humped nature. 

2 Environmentally coupled time-delayed 
system 

2.1 Environmental coupling 

Let us consider two first-order identical time-delayed 
systems coupled indirectly through a common environ- 
ment z. The mathematical model of the coupled system 
is given by 

x = -ax(t) + bif(x T ) + eifiiz, (la) 
y = -ay(t) + b 2 f(y T ) + ei&z, (lb) 

i = - —(Pix + ^22/). (lc) 

where a > 0, &i and b 2 are called the feedback rates 
for the x-system and ^-system, respectively Also, u T = 
u(t — t) (u = x,?/), where r G M + is a constant time 
delay e\ determines the coupling strength between en- 
vironment and the systems that controls the effect of 
environment on the systems. e 2 determines the coupling 
strength between the system and environment that con- 
trols the effect of individual systems on the environ- 
ment. Pi and fa determine the nature of coupling: for 
Pi = l,/^2 = —1 the systems are attractively coupled; 
on the other hand repulsive coupling is achieved for 
Pi = 1 and p 2 = 1. Finally k (> 0) determines the 
nature of the environment; in absence of both the x- 
and y- systems (i.e. e 2 =0), the environment decays 
towards the zero steady state and remains in that dor- 
mant state. 



2.2 System description and dynamics of the uncoupled 
system 

In this section we describe the time-delayed system pro- 
posed in Ref. [34 j, and also briefly discuss its important 
dynamical features. Ref. [34] proposed the following 
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Fig. 2 Phase plane plot in x-x(t — r) space for different r 
(6 = 1): (a) r = 1.40 (period-1), (b) r = 1.72 (period-2), 
(c) r = 1.94 (chaos), (d) r = 3.6 (double scroll hyperchaos). 
(other parameters are: a — 1, n — 2.2, m = 1, I = 10). 

first-order, nonlinear, retarded type delay differential 
equation with a single constant delay: 

x(t) = -ax(t) + bf(x T ), (2) 

where a > and b are system parameters. In [34] we 
consider both positive and negative values of 6, but in 
this paper we will consider only the case b > 0. Now, we 
define the following closed form mathematical function 
for the nonlinearity: 

f(x T ) — —nx T + ratanh(Zx r ), (3) 

where n, m and I are all positive system parameters and 
they are restricted by the following constraint: n < ml. 
It can be seen that the nonlinear function is constituted 
by the weighted superposition of two functions, namely, 
the linear proportionality function and the nonlinear 
tanh function. Further, f(x T ) is an odd-symmetric func- 
tion, i.e., f(-x T ) = -f(x T ). 

Fig |T] shows the nature of the nonlinearity produced 
by f(x T ) for different values of n, m and I. The non- 
linearity shows a hump in the first quadrant and the 
third quadrant. The condition n < ml ensures the two- 
humped nature of the nonlinear function. Also, Q has 
another distinct advantage - it provides a large num- 
ber of choices of n, m and I for which the two-humped 
nature will be preserved. 
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Fig. 3 The first eight Lyapunov exponents (A) with r; First 
two LEs become positive for r > 3.25 indicating hyperchaos. 
Other parameters are same as Fig. [5] 



The detailed stability and bifurcation analysis, ex- 
perimental implementation and results have been re- 
ported in [34]; there we have proved the existence of 
chaos and hyperchaos through the presence of strange 
attractor along with positive Lyapunov exponent [38] 
and higher values (> 3) of Kaplan- York dimension. 
Also we discussed the distinct features of this time- 
delay circuit over the existing piece-wise-linear nonlin- 
earity based time-delay circuits [39]. For better under- 
standing of the coupled case, let us discuss the impor- 
tant characteristics of the the system. Following pa- 
rameter values have been used in [34]: a = 1, n = 2.2, 
m = 1, I = 10. It has been shown that, for b = 1, if one 
varies r, for r > 1.102, the fixed point loses its stability 
through Hopf bifurcation. At r = 1.65, limit cycle of 
period- 1 becomes unstable and a period-2 (P2) cycle 
appears. Further period doubling occurs at r = 1.79 
(P2 to P4). Through a period doubling sequence, the 
system enters into the chaotic regime at r = 1.84. With 
further increase of r, at r = 2.60, the system shows 
the emergence of hyperchaos. The system shows a dou- 
ble scroll at r ~ 3.24. Phase plane representation in 
the representative x — x(t — r) plane for different r is 
shown in Fig (2] which shows the following character- 
istics: period-1 (r = 1.40), period-2 (r = 1.72), chaos 
(r = 1.94), and double scroll hyperchaos (r = 3.6). Fig- 
ure [3] shows the Lyapunov exponent spectrum of the 
system in the r parameter space. It is noteworthy that 
for a proper choice of 6, the system shows chaos and hy- 
perchaos even for a small time delay; e.g. for b = 1 one 
has chaos for r « 1.84 and hyperchaos for r > 2.60. 
This makes the circuit implementation of the system 
easier and also makes it superior for the possible appli- 
cations in communication system. 




Fig. 4 Experimental circuit diagram of the coupled system. 
R = 1 kft, C = 0.1 /zF, C 7 = 0.1 /xF. Buffers are designed 
with the unity gain non-inverting operational amplifiers. b x , y 
are amplifiers and ND represents nonlinear device (see Fig. 
[5}. A4x,y op amps are used as inverting adder and A5 op 
amp is used as inverting integrator. All the opamps are TL 
074. Rk — 1 kO, R& x ,y = R9x,y = 1 kO. In experiment the 
following condition is always met: Rio x = Rio y and Rn x = 

Rlly 



3 Experiment 

3.1 Electronic circuit realization 

We implement the coupled system given by in an 
analog electronic circuit. Figure [4] shows the represen- 
tative diagram of the experimental circuit. The pro- 
posed circuit consists three distinct parts, namely, the 
x-system, the y-system and the z-system or the envi- 
ronment The circuit of nonlinear device (ND) of each 
systems is given in Fig. [5] delay block is realized by 
using active all-pass filters (APF). To achieve the in- 
direct coupling via environment, we feed the outputs 
from the Ro — Co junctions of the x- and y- systems to 
the inverting terminal of the op- amp A5, which acts as 
an integrator, through buffers (for impedance match- 
ing) and resistors R\\ x and R\\ yi which determine the 
parameter 62- The voltage from the Ro — Co terminal 
of the ^/-system can be inverted for the in-phase cou- 
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pling (i.e. pi = 1,/^ = — 1) by the use of an unity 
gain inverter and connecting the points "D" and "E". 
Whether for inverse-phase coupling (Pi = 1,/^ — 1), 
terminals "D" and "F" will be connected. The output 
of the integrator A5 is fed into the inverting terminal 
of it through a buffer and a resistance Rk] Rk will de- 
termine the parameter n. Further, the output of the 
integrator is passed through an unity gain inverter and 
a buffer, and it is distributed in two ways: (i) it is di- 
rectly added to the x-system by the use of an inverting 
adder (A4x) through R\q Xi that determines the cou- 
pling strength ei, and (ii) it is passed through an in- 
verter (for Pi = l,/?2 = —1) or directly connected (for 
Pi = l,/?2 = 1) to inverting adder A4y through resis- 
tor Rioy that also determines ei; thus we always kept 
RiOx = Rioy and Rn x = Riiy 

Let, Vi(t) be the voltage drop across the capacitor 
Cq of the low-pass filter section of the x-system and 
that of the ^/-system be V 2 (t). Also, let the output of 
the integrator A5 be Vs(i), which represents the en- 
vironment. Thus the following equations represent the 
time evolution of the circuit: 



RqCq 



RqCq 



dVijt) 
dt 



dV 2 {t) 
dt 



^9x 
R~ 



lOx 



Rq 



-V 2 (t) + '-^b 2 f(V 2TD ) 



Rq 



dV 3 (t) 
dt 



RlOy 
1 



(4a) 



(4b) 



Rk 3 ^ ^ Rllx 



PiVi(t) 
1 



R 



lly 



foV 2 (t). (4c) 



Here, 61,62 = ^ is the gain of the amplifier A3i (i = 
x,y) (Fig(5|, Pi = 1, and P2 = =Fl, depending upon 
connection topology. f(V jTr> ) = f(Vj(t-T D ) = f(V jr ), 
(j = 1,2), is the nonlinear function representing the 
output of the Nonlinear Device (ND) of FigjHJ in terms 
of the input voltage Vj . Tjj is the time delay produced 
by the delay block. Also, we choose Rs Xj y — ^9x,y — 1 
kCl. In [34], we have reported that the nonlinearity of 
nonlinear device given by Fig. [5] has the following form: 



f(V jTD ) 



- -?*iv 



tanh( W fi|-Y 

V ttli V sat 



(5) 



Here P and w are the scaling factors that depend upon 
the non-ideal nature and asymmetry of the op amps. In 




Fig. 5 Nonlinear Device (ND) along with the gain block (bi). 
Ali-A3i (i = x,y) are op-amps (TL 074), R x = 1.26 kQ, 
R 2 = 19.29 kQ, R 3 = 52.81 kQ, R 4 = 6.73 kQ, R 5 = 10 kQ, 
Rq = 1 kf2. Inset shows the experimental oscilloscope trace 
of the nonlinearity produced by the ND. 




DELAY 



Fig. 6 Active first-order all-pass filter. R12 
C = 10 nF. "A" represents TL 074 opamp. 



Ri 



2.2 kfi, 



general, for smaller input voltages, P « 1 and w « 1. 
V S at is the saturation voltage of the opamp. The vari- 
able delay element is realized by a first order all-pass 
filter (APF) (Fig|6| [40 . The APF has the following 
transfer function: 



T(s) 



-CLl - 



UJ 



(6) 



with flat gain a\ = 1 (determined by Ri2i and ir^), 
and ujo = 1/CR is the frequency at which the phase 
shift is 7r/2. Since it has an almost linear phase re- 
sponse, thus each APF block contributes a delay of 
Td « RC. So n blocks produce a delay of = nRC 
(n = 1,2,...). By simply changing the resistance i?, 
one can vary the amount of delay; thus one can control 
the resolution of the delay line (the same technique of 
implementing delay line in the system and the coupling 
path has been used in [23) 120], which differs and is ac- 
tually advantageous over the conventional techniques 
used in j EB251E3 ) • 

Let us define the following dimensionless variables 
and parameters: t = R * c , r = J D ri , x = tt^, x t 
Vlr r- _ v 2 (t) _ V2r -u _ -.>vw _ 



v sat ■ 

R5i 

Rsi 



mi, w 



R2i — 

Rn 



v sat 



RqCq ' ^ V sat 

VsW Rsi . 
V sat ' R 4i - 

— R9i 

Rioi ' 



, z = 

— Rq 

Re>i 
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RkCii f = ^77c7' where i = x iV- To make the time " 
scale of the x, ^/-systems and the z- system equal we use 

Cj C . Now, we get k = f^, f = Thus, to 

make e\ = e 2l we have to use Rm = 2R\$i (since Rq — 

R 9i = 1 kfi). 

Now, the equations Q and ([5| can be reduced to the 
following dimensionless, coupled, first-order, retarded 
type nonlinear delay differential equations: 

dx 

— = -x(t) + bif(x T ) + ei^iz, (7a) 
^ = -y(t) + b 2 f(y T ) + e 1 f3 2 z, (7b) 

^ = - y + fay) • (7c) 

where, 

f(u T ) = — ni^ r + mi tanh(Ziu r ), (8) 
where u = x : y. 

It is worth noting that equations ([7]) (along with (J8|) 
are equivalent to equations ([!]) (along with (|3|) with 
a = 1, and appropriate choices of ni, mi and Zi. 

3.2 Experimental results 

The coupled system is designed in hardware level on 
a bread board using IC TL074 opamps (JFET quad 
opamps) with ±15 volt power supply. Capacitors and 
resistors have 5% tolerance . The resistance values used 
in the circuit for both the x- and ^-systems are: R\ — 
1.26 kft, R 2 = 19.29 kft, R 3 = 52.81 kft, R 4 = 6.73 
kO, = 10 kfi, = 1 kO. For the low pass sec- 
tion we used Rq = 1 kQ and Co = 0.1 fiF. Also for 
the z-system we choose Rk = 1 kQ and Cj = 0.1 /iF. 
The nonlinearity produced by the nonlinear device of 
both the x,y— systems are kept similar in nature and 
is shown in Fig. [5] (inset). The identical active all-pass 
filter stages of delay line (Fig.[6| have R\ 2 = -R13 = 2.2 
kf2, C = 10 nF and a variable resistance R. Here our 
main concern is to study the synchronization phenom- 
ena by varying the coupling strengths keeping the other 
system design-parameters same for the two systems. 

3.2.1 Real time waveform and phase plane plots 

To observe the real time waveform and phase plane 
plots we fix the feedback delay to the value r = 3.6 
(by choosing four stages of APF delay blocks with first 
three stages having R = 10 kO and the last one has 
R = 6 kfl). We choose Rj « 2.1 k^ to keep both the 
systems in the hyperchaotic regime. Figure [7] shows the 
hyperchaotic attractor of the x— system (upper panel) 




Fig. 7 The individual attractors of the uncoupled systems. 
Upper panel: x-system, Vi(t) — Vi(t — Td) space, lower panel: 
i/-system, Vz{t) — Vzit — To) space in the hyperchaotic regime. 
R7x,y =2.1 kO (for other parameter values see text). Oscil- 
loscope scale divisions: x -axis: 0.5 v/div, i/-axis: 0.5 v/div. 



and systems (lower panel) in the uncoupled condi- 
tions. In the following cases, we consider the coupling 
resistors for which e\ = e 2 . 

(i) For fa = 1 and fa = —1: This condition is satisfied 
by connecting the points "A" to "B" , and "D" to "E" , 
in the Fig. [4j Figure [8] shows three distinct cases: (a) 
The x- and ^/-systems evaluate independently in time 
for large values of coupling resistance. For example, 
when we keep R\q x — Rio y ~ 9.96 kO, R\\ x = R\\ y 
20 kfi, both the systems are unsynchronized. Fig.[8^al) 
shows the time series of the x-system (yellow) and the 
^/-system (blue), and the red trace in the lower portion 
of this plot shows the difference (V\{t) — V 2 (t)) (call 
it the "error") of the two systems captured by a digi- 
tal storage oscilloscope (DSO) (Tektronix TDS2002B, 
60 MHz, 1 GS/s). It can be noticed that the ampli- 
tude of the error signal is of the same order as the 
original signals V\(t) and V 2 (t)) indicating unsynchro- 
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Fig. 8 fi\ = 1 and p2 — — 1: Experimental time series plot 
of the x-system Vi (t) (yellow) and the ^-system V2 (t) (blue) 
in the hyperchaotic Regime, lower trace in red represents the 
error signal (Vi(t) — V^it)): (al) unsynchronized state, (bl) 
in-phase synchronization (cl) complete synchronization. The 
corresponding phase plane plots are shown in (a2), (b2), and 
(c2), respectively. (For parameter values see text). For (al), 
(bl), and (cl): x-axis: 25/xs/div, y-axis: 1.25 v/div. For (a2), 
(b2), and (c2): x-axis: 1 v/div, y-axis: 1 v/div. (Color figure 
online). 






V 

\ 



\ 



Fig. 9 /3i = 1 and ^2 = 1: Experimental time series plot 
of the x-system Vi (t) (yellow) and the i/-system V2 (t) (blue) 
in the hyperchaotic Regime, lower trace in red represents the 
error signal (Vi(t) + V^(t)): (al) unsynchronized state, (bl) 
inverse-phase synchronization (cl) anti-synchronization. The 
corresponding phase plane plots are shown in (a2), (b2), and 
(c2), respectively. (For parameter values see text). For (al), 
(bl), and (cl): x-axis: 25 /xs/div, y-axis: 1.25 v/div. For (a2), 
(b2), and (c2): x-axis: 1 v/div, ^/-axis: 1 v/div. (Color figure 
online). 



nized states. Fig. |8|a2) shows the corresponding phase- 
plane plot in Vi (t) — V2 (t) space that confirms the asyn- 
chronous condition, (b) If we decrease the coupling re- 
sistors (i.e., increase the coupling strength e), the cou- 
pled systems show in-phase synchronization. This is 
shown for R\q x = Rio y ~ 1.57 kO, R\\ x = Rn y ~ 3.2 
kQ in Fig.|Hbl)(b2); the time series m Fig.^bl) shows 
that the waveforms are in phase but their amplitude 
levels are still uncorrelated in time scale. Also the am- 
plitude of error signal is much reduced here. Fig. (8^b2) 
shows the corresponding phase-plane diagram where we 
can see that the system dynamics now wanders around 
the 45° diagonal line . These two indicate the in-phase 
synchronization of the x- and the ^/-systems, (c) Fur- 
ther reduction in the coupling resistors results in com- 
plete synchronization. Fig. [8] (cl)(c2) shows this for 
RiQx = R\<dy ~ 0.56 kfl and R\\ x = R\\ y ~ 1.2 kfl. 
Fig.[8jcl) and (c2) shows the real time and phase-plane 
plots of the systems, respectively. The error line in lower 
portion of Fig. (8^cl) indicates that the two waveforms 
are equal in phase and amplitude. The phase-plane plot 



in Fig.[8jc2) showing 45° inclination with both the axes 
confirms complete synchronization. 

(ii) For P\ = fa = 1: This condition is achieved 
when one connects the points "A" to "C" and "D" to 
"F" , in Fig. [4] In this case also we get three distinct 
situations: (a) The x- and y- systems evaluate inde- 
pendently in time for larger coupling resistance. Figure 
|9jal)(a2) show this for Riq x — Rio y ~ 9.93 kfi, Rn x = 
Riiy — 20 k£l. Fig.[9^al) shows the time-series of the x- 
and systems, and that in the lower portion of the plot 
shows the sum of the two waveforms (Vi(t) + V2(t)). 
Fig. |9ja2) shows the corresponding phase-plane plot. 
The sum and the phase-plane plot confirm that the 
systems are not synchronized, (b) Lowering of cou- 
pling resistances results in inverse-phase synchroniza- 
tion. This is shown in Fig.^bl) (b2) for Rio x = Rio y = 
2 kft, R llx = R lly « 4.1 kft. Fig. [9pol) shows the 
time-series and the sum (Vi(t) + V<2(£)), and Fig. |9jb2) 
shows the corresponding phase-plane plot. One can see 
that Vi(t) and V^(t) are in the phase- inverted mode 
but their amplitudes do not correlate, (c) At further 
lower coupling resistance values the systems show anti- 
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synchronization; Fig.[9|cl) and (c2) shows the real time 
and phase-plane plots, respectively, for R\q x = Rio y ~ 
0.3 k£7, Rii x = Riiy « 0.7 kfi. Here one can see that 
the two systems have ir phase shift and also their ampli- 
tude levels are same. The phase-plane plot shows that 
the system dynamics lives in a diagonal line making 
an angle 135° with both the axes. This confirms the 
occurrence of anti-synchronization. 

3.2.2 Generalized autocorrelation function and 
correlation of probability of recurrence 

It is seen from numerical simulations (Fig. |5| and ex- 
perimental results (Fig. [7]) that the attractor of the 
system under study is not phase-coherent. As there ex- 
ists no general technique of finding phase of a phase- 
incoherent hyperchaotic attractor, we use two dynami- 
cal measures of phase synchronization (PS) proposed in 
Ref. [35j[36], namely, generalized autocorrelation func- 
tion (P(t)) and correlation of probability of recurrence 
(CPR). These two measures have effectively been used 
in the context of chaotic phase synchronization of cou- 
pled time-delayed systems (e.g. [25], [27], [41], etc.). 
Here, we compute P(t) and CPR of the coupled sys- 
tem experimentally using the time-series data acquired 
from the experimental circuit. 

The generalized autocorrelation function (P(t)) is 
defined as 35,36 



Pit) 



N 



N-t 



WXi - x, 



i+t\ 



(9) 



where is the Heaviside function, is the i — th 
data point in the X variable, N is the total number 
of data points, e t is a preassigned threshold value, and 
||.|| represents the Euclidean norm. Let, P\{t) repre- 
sents the generalized autocorrelation function of the x- 
system and be that of the ^/-system. We compute 

P\{t) and P2(t) from the experimental time-series data 
(N = 2400) acquired using DSO (Tektronix TDS2002B, 
60 MHz, 1 GS/s). For both the cases we choose the 
threshold value e t = 0.01. 

1,02 



Fig. 10 shows P(t)s for fii 



-1 (i.e. in- 

phase coupling case), and Fig. 11 shows the same for 
01 = 1, 02 = 1 (i-e. inverse-phase coupling case). Upper- 
panel of Fig. [To] represents the unsynchronized state for 
01 = 1,^2 = -1 at R 10XiV « 9.96 kfi, R 11XtV = 20 
kfl (same parameter values as of Fig{8^al)); similarly 
Fig. 11 (upper-panel) shows the same for /3i = 1, 02 = 1 
at Riox, y ~ 9.93 kfi, Ru x ,y = 20 kQ (same parameter 
values as of Figj9jal)). From both the figures we can see 
that peaks of P\(t) does not match with that of 
in the t-axis, confirming the fact that the phases of 
the two oscillators are not synchronized. Lower panels 



of Fig. [10] (R 10XtV w 1.57 kfi, 

Fig. 11 (R\0x,y = 2 kfi, Rll x ,y ' 



Riix,y ~ 3.2 kO) and 
^4.1 kO) show that the 
dominant peaks of P\(t) and P2(t) matches exactly in 
the t-axis confirming the phase synchronization of the 
coupled oscillators. Also, almost equal amplitude of the 
peaks in the phase synchronized cases shows the good 
quality of PS in the coupled systems. 

A quantitative measure of P(t) is the correlation 
of probability of recurrence (CPR) that was defined in 
Ref. [35][36] as 



CPR = 



(Pi(t)P 2 (t)) 



<7i<72 



(10) 



A, 2 present that the mean value has been subtracted, 
and <t 1?2 are the standard deviations of the P\{t) and 
P2W, respectively. In the phase synchronized case, gen- 
erally, CPR^ 1, and for the unsynchronized cases its 
value is appreciably smaller than 1. Using the definition 



of ( 10 ), we compute CPR from experimental time-series 



data both for the in-phase and inverse-phase cases. For 
in-phase synchronized case we have CPR=0.997 (pa- 



rameters are same as of lower panel of Fig. 10), and for 



inverse-phase synchronized case we have CPR=0.985 



(parameters are same as of lower panel of Fig. 11), 
which confirm the occurrence of phase synchronization 
in the coupled system. 

3.2.3 Concept of localized sets 



We use another dynamical measure for qualitative con- 
firmation of phase synchronization called the concept of 
localized sets (CLS) proposed in [37 . It was shown in 
[37] that the CLS technique is extremely useful to detect 
phase synchronization even when no proper measure of 
phase is possible. The idea of CLS in coupled oscilla- 
tors is based on the fact that if one identifies a certain 
event in the first oscillator and then track the second 
oscillator at that particular event, a set D will be ob- 
tained for the second oscillator; if that set spreads over 
the whole attractor space of the second oscillator then 
one may say that there is no phase correlation between 
the two coupled oscillators. On the contrary, if the set 
D becomes localized to a certain zone of the attractor 
space then one can say that the coupled oscillators are 
phase synchronized. 

From the experimental circuit we simultaneously ac- 
quire the time-series data of x- and ^/-systems (with 
2400 data points for each). Next, we define the event 
Vi(t) — —0.5 for the x— system and track the values of 
Vi(fy from the time-series whenever that event is met. 
The obtained set of data values of T^(t), and the corre- 
sponding V\ and V2 are plotted in V\ — V2 space. Fig- 
ure. [12] shows the localized sets for the in-phase and 
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t 

Fig. 10 ft = l,/?2 = — 1: Generalized autocorrelation func- 
tion (-P(t)) computed from the experimental time-series of 
x-system (Pi, in red ) and ^-system (P2, in blue) in the un- 
synchronized state (upper panel) and in-phase synchronized 
state (lower panel). Notice the matching (difference) of the 
peaks of Pi and P2 in the in-phase synchronized (unsynchro- 
nized) case (for parameter values see text). (Color figure on- 
line). 




100 200 300 400 500 



Fig. 11 Pi = 1, P2 = 1: Generalized autocorrelation func- 
tion (P(t)) computed from the experimental time-series of 
^-system (Pi, in red) and y-system (P2, in blue) in the un- 
synchronized state (upper panel) and inverse-phase synchro- 
nized state (lower panel). Notice the matching (difference) 
of the peaks of Pi and P2 in the inverse-phase synchronized 
(unsynchronized) case (for parameter values see text). (Color 
figure online). 



inverse-phase cases at unsynchronized, phase synchro- 
nized and complete (anti-) synchronized cases. It can 
be seen that for the unsynchronized cases ( Fig [l2] (a) 
for ft = 1,02 = -1 and Fig. [12] (d) ft = l,ft = 1), 
the set D (represented by black points) spreads over the 
whole attractor space of V2 (t) (resistor values are same 
as used in Fig.[8](al) (ft = l,ft = -1) and Fig. [9] (al) 
(ft = 1, ft = 1), respectively). With further decrease in 
Ri0x,y and Riix, y (that is increase in coupling strength) 
the set D becomes localized in V2(£)-space indicating 
the occurrence of phase synchronization; Fig.[l2](b) and 
(e) show this case for in-phase and inverse-phase cases, 
respectively (resistor values are same as used in Fig. [8] 
(bl) (ft = 1,02 = -1) and Fig.|9](bl) (ft = l,ft = 1), 
respectively). Here, we can see that the black points 
representing the set D are localized to a small zone 
of V^(t). Another interesting observation can be made 



from Fig. [l2j that is, for in-phase synchronized state 
(Fig. [l2] (b)), the localized set D is situated around 
Vzii) = —0.5, which is equal to V\{t) = —0.5, whereas 



for inverse-phase synchronized state (Fig. 12 (e)), the 
localized set D is situated around V2(t) — 0.5, which is 
opposite of Vi(t) = —0.5 (indicating a tt phase shift be- 
tween V\(t) and V2(t)). Finally, with further decrease in 
coupling resistance (i.e. increase in coupling strength) 
we observe that the set D becomes localized to a very 
narrow range of V2 (t) that indicates complete synchro- 
nization (Fig. [l2] (c)) and anti-synchronization (Fig. [l2] 

(f))- 



4 Linear stability analysis 

In this section we explore the linear stability of syn- 
chronized states of the coupled system of 0. In [29] 
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(e) 





Fig. 12 Concept of localized sets computed from the experi- 
mental time-series: Black dots represent the set D for the par- 
ticular event V\(t) = —0.5, the points are plotted along with 
Vi(t) and V 2 (t) in the Vi(t) - V 2 (t) space. p 1 = I, f3 2 = -1: 
(a) unsynchronized (b) in-phase synchronized (c) complete 
synchronized states. f3\ = 1,(32 — 1: (d) unsynchronized (e) 
inverse-phase synchronized (f) anti-synchronized states. Note 
the localization of the set D in the narrow region in V2 (t) axis 
for the phase synchronized (b,e) and complete (anti-) syn- 
chronized cases (c,f). For parameter values see text. (Color 
figure online). 



stability of the synchronized states of the environmen- 
tally coupled flows with no delay has been derived with 
some broad approximations; It was shown there no ex- 
act analysis is possible; nevertheless, the authors ar- 
rived at a condition that predicts the stable complete 
(anti-) synchronized zone in parameter space. In the 
present case the scenario is more complex owing to the 
presence of delay. Let us start by considering -0, and 
<fi be the deviations from the synchronized states of the 
system variables in the Eq. 0. Then linearizing the 
system along with these deviations leads to 



-atp(t) + bif{x T )tp T + ei/?i<£, 
-a0(t) + b 2 f'(y T )6 T + £l /? 2 0, 



€2 



<t> = -K4--Z{M + h0)- 



(lla) 
(lib) 
(11c) 



where, u T = u(t — r), u = ip, 0. 

It is not possible to carry out an exact analysis of 



(11). To make the analysis possible we impose some 



constraints to it and consider this one as a special case. 
Let us consider the complete synchronized state of the 
systems, i.e, x = y and hence, x T = y T . Also we define 
a new variable relating ip and in the following way: 



(12) 



With this, Eq. ([TTJ can be reduced to following equa- 
tions (with bi = bi = b): 



X = -a X + bf(x T )xr + + PDh 
£2 



<j) = — K<j) 



-x- 



(13a) 
(13b) 



Equation (13) also can not be solved in closed form. We 



can make further approach by considering f'(x T ) = 5' ', 
where 8' is a constant; this approximation was also used 
and justified in [29 . Define S = b6', and $\ + 0% =2, 
from ( 13 ) we get 



X 



(j) = —KC)) 



^X + ^Xt + 2ei</>, 
£2 



-x- 



The characteristic equation of ( 14 ) is given by 
Se~ Xr 



det 



A + a 



C2 

2 



"26! 
A + K 



0, 



which gives on evaluation 

A 2 + (a + K - Se- XT )\ - 5Ke~ XT 



+ (an + eie 2 ) = 0. 



(14a) 
(14b) 



(15) 



(16) 



The eigenvalue of the characteristic equation ( 16 ) may 



be real or imaginary. Let us consider A = \i ± iv. The 
synchronization will become just oscillatory if the eigen- 
value be a purely imaginary one. Thus for the limiting 
case, we consider \i = 0, and A = ±.iv. Substitution of 
this in the above and a comparison between the real 
and the imaginary parts yields 

2 



bS ' v sin vr — bS' k cos vt 
-bS'is cos vr + bS'ft sin vr 



v" - {an + 6ie 2 ), 
-(a + k)v. 



(17) 
(18) 



Squaring and adding Eq. (17) and (18), we get 



^ 4 + (a 2 + ^ - 2exe 2 - 6V V 



+ (a^ + eie 2 ) 2 - 6 W = 0. 



(19) 



2eie 2 -6 2 (5 /2 ) = O and (a^ + eie 2 ) 2 



Consider (a 2 + , 

b 2 5' 2 n 2 = A. Thus Eq. ([l9]) reduces to 



<9^ 2 +yl = 0. 



Eq. (20) has the following solution: 



-O ± V<9 2 - AA 



(20) 



(21) 
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From the above equation it may be stated that v 2 must 
be real and positive, otherwise there is no purely imagi- 



nary roots of Eq. ( 16 ), and for this one requires Q 2 > AA 
along with v 2 > 0, that results 



6162 < \k(S — a)\. 



(22) 



This condition leads to the oscillatory solution of the 
complete synchronization (CS) indicating the loss of 
CS. Thus, the threshold condition of CS is 



(23) 



Here, e\ cr and €2 cr are only the lower limits of the 
coupling strengths where a stable complete synchro- 
nized state can be achieved. If we increase the coupling 
strength beyond the critical values, then the synchro- 
nized state will be prevailed. In the next section we nu- 
merically confirm this result with proper choice of the 



effective value of 5. Also note that (23) is equally valid 
for the anti-synchronized (AS) state, because in that 
case consideration of y(t) = —x(t) and f3\ = = 1 
does not alter the form of Jill). 



5 Numerical Simulation 

The system equation is simulated numerically us- 
ing Runge-Kutta algorithm with step size h = 0.01. 
Following initial functions have been used for all the 
numerical simulations: for the x-system (j) x (t) = 1, for 
the ^/-system <j> y (t) = 0.9, and for the environment, i.e., 
z-system, z(0) = 0.85. Also, the following system de- 
sign parameters are chosen throughout the numerical 
simulations: a = 1, n = 2.2, m = 1, I = 10, k = 1; also 
we choose b\ = 62 = 1, and r = 3.6 to ensure that both 
the systems be in the hyperchaotic region [34]. Unless 
otherwise stated, we use e\ = €2 = e. In all the com- 
putations, and real time and phase plane diagrams, a 
large number of iterations have been excluded to allow 
the system to settle to the steady state. 



5.1 Lyapunov exponent spectrum 

One can detect the onset of in-phase (inverse-phase) 
synchronization and complete (anti-) synchronization 
directly from the Lyapunov exponent (LE) spectrum of 
the coupled system [42]. We compute the LE spectrum 
in e-space (ei = 62 = e) directly from 0. Among a 
large no of LEs we track the behavior of the first four 
LEs that are sufficient to detect the occurrence of phase 
or complete (anti-) synchronization with the variation 
of e [29 . For ft = l,/3 2 = -1, Fig. [l3] (a) shows the 
first five LEs in e-space; it can be observed that with 
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Fig. 13 Lyapunov exponent spectrum of the coupled system, 
(a) £1 = l,/3 2 = -1. (b) £1 = l,/3 2 = 1. Transitions of LE4 
from positive to negative value through zero indicates phase 
synchronization, similar transition of LE3 indicates complete 
(anti-) synchronization. 



increase in e, the fourth largest LE (LE4) crosses the 
zero value to become negative at e ~ 1.1 indicating 
the transition from unsynchronized state to in-phase 
synchronized state. With further increase in e, we ob- 
serve that LE3 also becomes negative (from a positive 
value) at e ~ 1.48, which indicates the transition from 
in-phase synchronization to complete synchronization. 



Similar case is shown in Fig. 13 (b) for /3i = l,/?2 = 1; 
here we can see that LE4 becomes negative at e ~ 1.2 
indicating inverse-phase synchronization and transition 
of LE3 from a positive value to a negative value at 
e w 1.56 indicates anti-synchronization of the coupled 
systems. 



5.2 Time series and phase-plane plots 

Guided by the above result, we choose three different 
coupling strengths, e (ei = €2 = e), and plot the time 
series and phase-plane plots for the illustrative exam- 
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pies of unsynchronized, phase synchronized and com- 
plete (anti-) synchronized states. Figure. [l4| shows this 
for ft = l,ft = — 1- Upper panel of Fig. [14| shows the 
unsynchronized case at e = 0.2; corresponding phase- 
plane plot indicates no correlation between the x-system 



and ^-system. Middle panel of Fig. 14 shows the in- 
phase synchronized state at e = 1.4; we can see from 
the corresponding phase-plane plot that although the 
phases of two systems are synchronized but their am- 



plitudes still differ. Finally, lower panel of Fig. 14 shows 
the case of complete synchronization at e = 1.6; here 
both the systems become identical to each other. 
Figure. [15] shows this scenario for ft = ft = 1. At 
e = 0.2 we can see the unsynchronized case (upper 
panel); at e = 1.4 the system is inverse-phase synchro- 
nized (middle panel), and finally at e = 1.6 we observe 
anti- synchronization (lower panel), where x(t) and y(t) 
have a 7r phase shift and their amplitude levels are sim- 
ilar. 



5.3 Generalized autocorrelation function and CPR 



As discussed in Sect 3.2.2 we compute P(i) and CPR 
from numerical simulations. We use N = 5000 and 
plot Pipit) for the unsynchronized and phase synchro- 
nized cases. Figure fl6| (a) and Fig|l7| (a) shows this at 
ei = e 2 = 0.2 for (ft = 1, ft = -1) and (ft = ft = 1), 
respectively. The unmatched dominant peaks of Pip(t) 
indicates the lack of phase synchronization in the cou- 
pled systems at this small value of coupling strength. 
We plot the same at 61=62 = 1.4 (Figure 16 (b) 
and Fig [l7] (b)). Here the perfect matching of domi- 
nant peaks of Pipit) indicates the occurrence of PS; at 
this point we computed CPR, which is equal to 0.99 for 
Pi = l,ft = -1, and 0.987 for ft = ft = 1 that also 
confirms the occurrence of PS in the coupled system. 



5.4 Concept of localized set 

Further, we use the technique of concept of localized 
set to numerically confirm PS in the coupled system. 
Here, we define the following event in the y- system: 
y(t) = —0.2 and y(t — r) < —0.2; whenever this event 
occurs in time, we track the values of the x-system and 
get a set of values of x(t). Then we plot this set in the 
x(t) — x(t — r) space. The similar process is repeated for 
the ^/-system with the following event in the x-system: 
x(t) = —0.2 and x(t — r)< —0.2. Figure. 18 shows this 
for ft = l,ft = —1 and Fig |l9] represent the case for 
Pi — P2 — 1 for various coupling strengths. Figure. 18 
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Fig. 14 Pi = 1 and ft = — 1: Numerically obtained time- 
series of x and y, and corresponding phase-plane plots in x — y 
space. Upper panel: unsynchronized state (e = 0.2), middle 
panel: in-phase synchronized state (e = 1.4), lower panel: 
complete synchronized state (e = 1.6). (Color figure online). 



(a,d) and Fig. n9\ (a,d) show the spreading of data set 
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Fig. 15 f3i — 1 and ^2 = 1: Numerically obtained time-series 
of x and and corresponding phase-plane plots in x — y space. 
Upper panel: unsynchronized state (e = 0.2), middle panel: 
inverse-phase synchronized state (e = 1.4), lower panel: anti- 
synchronized state (e = 1.6). (Color figure online). 
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Fig. 16 (3i = 1,(32 — —1: Numerically obtained Generalized 
autocorrelation function, P(t), of x-system (Pi, in red) and 
i/-system (P2, blue) in the unsynchronized state (e = 0.2) (up- 
per panel) and in-phase synchronized state (e = 1.4) (lower 
panel). Notice the matching (difference) of the peaks of Pi 
and P2 in the in-phase synchronized (unsynchronized) case 
(for other parameter values see text). 



over the whole phase space for e\ = 62 = 0.2 indicat- 
ing phase- incoherent behavior of the coupled systems. 
At 6\ = 62 = 1.4 the data set forms a localized set in 
the phase space indicating the occurrence of phase syn- 



chronization (Fig. |18| (b,e) and Fig. 19 (b,e)). Finally, 
at ei = €2 = 1.6 we can see that Fig. [l8] (c,f) shows 
that the localized set becomes identical with the defined 
event itself for j3\ = l,/?2 = — 1, indicating a complete 
phase synchronization. Whereas, for /3x = ^ = 1, at 
ei = 62 = 1.6 Fig. [l9| (c,f) shows that the localized set 
is situated in a narrow region that is just opposite to the 
defined event, indicating anti-phase synchronization. 

5.5 Stability of synchronization in parameter space 
To support the theoretical result obtained analytically 



in (23), numerically we study the synchronization sce- 
nario in the e\ — 62 space. We use all the above men- 
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Fig. 17 Pi = l,/?2 = 1: Numerically obtained Generalized 
autocorrelation function, P(t), of x-system (Pi, in red) and 
i/-system (P2, in blue) in the unsynchronized state (e = 0.2) 
(upper panel) and inverse-phase synchronized state (e = 1.4) 
(lower panel). Notice the matching (difference) of the peaks 
of Pi and P2 in the inverse-phase synchronized (unsynchro- 
nized) case (for other parameter values see text). 








Fig. 18 Pi = 1 and P2 — — 1: Numerically computed Con- 
cept of localized sets: (a), (b), and (c): plot of x(t) — x(t — r) 
along with the set D (black dots) obtained by defining the fol- 
lowing event in the i/-system: y(t) = —0.2 and y{t — r) < —0.2 
(blue line in (d)). (d), (e), and (f): plot of y(t) — y(t — r) along 
with the set D (black dots) obtained by defining the following 
event in the x-system: x(t) — —0.2 and x(t — r) < —0.2 (blue 
line in (a)), (a) and (d) represents unsynchronized state with 
e — 0.2 (note that the black dots are scattered all around 
the attractors representing that the attractors are incoher- 
ent), (b) and (e) represents phase synchronized states with 
e = 1.4. (c) and (f) represents complete synchronized case 
with e = 1.6. Note that in the last two cases black dots form 
a localized set that lives in a narrow region of the attrac- 
tor. (Color figure online). 



tioned measures to detect the complete-synchronization 
(CS) and anti-synchronization (AS). The threshold value 
of ei and 62 for which CS and AS occurs are plotted 



in Fig 20 (a) and (b), respectively. The solid line rep- 
resents the theoretical prediction and points represent 
the numerically obtained threshold values for which CS 
or AS just start. In theoretical curve, for CS, we use 
5' = 3.20 (also a = 1, 6=1, and n — 1); this value 
of 5' is obtained for x T = ±0.0825 (remember that 
5' = f'(x T )). Similarly for AS we use 8' = 3.15, that 
implies x T ~ ±0.0832. In both the cases the values of 
x r , for which theoretical and numerical results match, 
lie well within the actual phase space where the system 
dynamics live. Thus the stability analysis agrees with 
the numerical results with some effective choice of 5' . 
Since the stability analysis relies on some broad approx- 
imations it requires further refinements and scopes are 
there to progress in that direction. 



6 Conclusion 

In this paper we have reported different synchronization 
scenario of hyperchaotic time-delayed systems coupled 
indirectly through a common environment. The system 
we have chosen is a first-order, nonlinear time-delayed 
system that posses a closed form mathematical function 
for the nonlinearity, shows hyperchaos even at a mod- 
erate or small time delay, and is convenient for the elec- 
tronic circuit design. For the first time we have explored 
the experimental aspects of environmentally coupled 
time-delay systems. We have confirmed the occurrence 
of phase synchronization and complete (anti-) synchro- 
nization in the experimental circuit by using recurrence 
analysis and the concept of localized sets computed di- 
rectly from the experimental time-series data. We have 
shown that with the proper choice of coupling param- 
eters (i.e. Pi = l,p2 = —1), increase in the coupling 
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Fig. 19 (3i — 1 and p 2 — \\ Numerically computed Concept 
of localized sets: (a), (b), and (c): plot of x(t) — x(t — r) along 
with the set D (black dots) obtained by denning the following 
event in the ^-system: y(t) = —0.2 and y(t — r) < —0.2 
(blue line in (d)). (d), (e), and (f): plot of y(t) — y(t — r) 
along with the set D (black dots) obtained by defining the 
following event in the x-system: x(t) = —0.2 and x(t — r) < 
—0.2 (blue line in (a)), (a) and (d) represents unsynchronized 
state with e = 0.2 (note that the black dots are scattered 
all around the attractors representing that the attractors are 
incoherent), (b) and (e) represents phase synchronized states 
with e = 1.4. (c) and (f) represents anti-synchronized case 
with e = 1.6. (Color figure online). 

strength results in a complete synchronized state from 
unsynchronized states via phase synchronized states. 
Also, for Pi = p 2 = 1 5 with the increasing coupling 
strength we have observed a transition from unsynchro- 
nized states to the anti-synchronized state via inverse- 
phase synchronized states. To corroborate the experi- 
mental results we have presented a linear stability anal- 
ysis of the complete (anti-) synchronized state, and 
also we perform detailed numerical simulations. Numer- 
ically we computed the Lyapunov exponent spectrum 
of the coupled system from where we have identified 
the zone of phase synchronization and complete (anti-) 
synchronization in the parameter space. We have used 
the recurrence analysis and the concept of localized set 
to numerically reconfirm the occurrence of phase syn- 
chronization in the system. Finally, using all the mea- 
sures we have identified the parameter zone in the two- 
parameter space for the stable complete and anti syn- 
chronization, which agree well with the theoretical re- 
sults. 
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Fig. 20 Stability curve in e± — e 2 parameter space (a) Pi = 1 
and (3 2 = -1, (b) (3i = 1 and (3 2 = 1. 

The present study can be extended to explore the 
synchronization scenario of other environmentally cou- 
pled time-delay systems, e.g. Mackey- Glass system, Ikeda 
system, etc. One of the main features of the environ- 
mental coupling scheme is that here the system dynam- 
ics in the synchronized states and unsynchronized states 
have almost the same structure in phase-space [29 ; this 
particular feature makes this coupling scheme advanta- 
geous over other time-delay coupling schemes. From the 
academic interest, as the environmental coupling is very 
much relevant to biological systems, and delay is inher- 
ent in biological processes, thus the present study will 
be important for the study of synchronizations in the 
biological systems with delay. Apart from the academic 
interest, the present study is important from applica- 
tion point of view also; for example, since with environ- 
mental coupling one can make a transition from com- 
plete synchronization to anti-synchronization by sim- 
ply changing the sign of /?2, thus this will be useful 
in easy implementation of hyperchaotic binary-phase- 
shift-keying (BP SK) -based digital communication sys- 
tems. 
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